Multiple 6-jets reveal natural SUSY and the 125 GeV Higgs 
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We demonstrate that in natural supersymmetry (SUSY), h — > bb can be a very promising channel 
to discover the SM-like Higgs resonance at the Large Hadron Collider (LHC), although in general 
contexts it is thought to be challenging because of its huge QCD background. In this scenario, 
the SM-like Higgs boson is mainly produced via cascade decays initiated by pair-produced stop or 
sbottom squarks. The good sensitivity to h — > bb owes a great deal to the application of multiple 
(> 4) b-jet tagging in removing the QCD background, and color-flow variables for reconstructing 
the Higgs resonance. We show in two benchmark points that a SM-like Higgs resonance can be 
discovered at the 14 TeV LHC (with a signal-to-background ratio as high as 0.3 — 0.4), with < 40 
fb _1 of data. Potentially, this strategy can be also applied to non-SUSY theories with natural EW 
symmetry breaking and cascade decays of top partners for the SM-like Higgs search, such as little 
Higgs, composite Higgs, and Randall-Sundrum models. 



Deciphering the Higgs mechanism is one of the top 
priorities of the ATLAS and CMS experiments at the 
Large Hadron Collider (LHC). Recently, both the CMS 
and the ATLAS collaborations announced discovery of a 
Higgs-like resonance, based on a combined analysis for 
the Standard Model (SM) Higgs searches via h — > 77 
and h — > ZZ* — > Al, with the reconstructed invariant 
mass ~ 125 — 126 GeV [l], Q- Furthering the identifica- 
tion of this new particle greatly enhances the necessity 
and emergency of information from complementary Higgs 
search channels like h — > bb. 

Unlike h — > 77 and h — > ZZ* -> 41, h -> bb can help 
test the mass generation mechanism of the SM fermions 
directly. However, in spite of its relatively good sensi- 
tivity at the Tevatron, the h — > bb search at the LHC 
is challenging because of its huge QCD background, ex- 
cept in the highly boosted regime where jet kinematics 
allows for the successful application of jet-substructure 
tools 0, 0]. Physics beyond the SM can modify the 
collider phenomenology of the Higgs search modes if 
there are new Higgs production mechanisms available. 
With new strategies designed, the sensitivity of measur- 
ing the challenging decay modes could be sizably im- 
proved. In this letter we will demonstrate that in a series 
of well-motivated scenarios like natural supersymmetry, 
the Higgs discovery via h — > bb can be greatly assisted by 
a set of dedicated strategies. 

Supersymmetry is the prime candidate theory for 
physics beyond the SM. It provides a natural solution to 
the hierarchy problem. In the scenario of natural SUSY, 
the scale for the parameters closely related to EW sym- 
metry breaking (EWSB) is not sizably higher than the 
EW scale, leading to a small fine-tuning in the EWSB. 
Compared to general SUSY scenarios, natural SUSY is 
theoretically more predictive and experimentally more 
accessible. It predicts the existence of light superpar- 
ticles [H|: two stop and one left-handed sbottom squarks, 
with their masses < 700 GeV; two Higgsino-like neutrali- 
nos and one Higgsino-like chargino, with their masses 
< 350 GeV; and a gluino, with its mass < 1500 GeV. 



With conserved R-parity, stop and sbottom squarks 
are pair-produced at colliders. An interesting feature is 
that they often decay into a SM-like Higgs boson, in asso- 
ciation with a b quark (for discussions on SUSY-assistcd 
Higgs production in various contexts, see 0-0] )■ If the 
Higgs boson decays into a pair of b quarks, the final state 
typically contains at least four b quarks, independent of 
the details of the decay chains. There are two main mech- 
anisms for the Higgs production via the stop and sbot- 
tom cascade decays. In the first case, the Higgs boson 
is produced via neutralino or chargino decays where the 
Yukawa couplings of the SM fermions are not directly in- 
volved. In the second case, the Higgs boson is directly 
produced via the decay of the heavier stop squark where 
the top Yukawa coupling gets involved. Because these 
processes are initiated by light squarks, the produced 
Higgs bosons tend to be less boosted, compared to the 
case discussed inl8(] (for discussions on less boosted Higgs 
boson, also see [7J). Two topologies of these mechanisms 
are illustrated in Fig. [T] Given that the QCD events 
containing multiple 6-jets are few, this feature provides a 
new strategy of improving the sensitivity of the h — > bb 
resonance search at the LHC. 
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FIG. 1: The primary Higgs-producing SUSY cascades in the 
first (left) and second (right) benchmark scenarios. 

To design the minimal strategies that ensure optimal 
coverage of the full space of models, we perform a col- 
lider analysis in two benchmark scenarios which repre- 
sent these two mechanisms, respectively. Their parti- 
cle mass spectrum and decay branching ratios are calcu- 
lated using SUSY-HIT @, in the context of the MSSM 
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with mh = 125 GeV set as an input (see Table U). Here 
new contributions to rrih have been implicitly assumed 
which could be either from new F-terms or from new 
D-terms (for recent discussions, e.g., see [loj|). These 
benchmarks are not fully realistic, but simplified mod- 
els in essence [ll|. We tolerate a [i parameter that is 
slightly higher than is usually desired in natural SUSY 
in our second benchmark point, but this is not neces- 
sary for the application of our methods. In any case, we 
do not adhere to a precise measurement of naturalness. 
Though searching for experimentally allowed points in 
the parameter space of natural SUSY is not the main 
motivation of this letter, we would like to point out that 
these benchmark points arc still pretty safe, after com- 
paring them with the publicly available searches, such as 



the CMS 6-enriched razor- variable search |L2j, the AT- 
LAS light sbottom search 0, and the CMS 6-jets +$ T 
search 
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Benchmarks 


I (GeV or %) 


II (GeV or %) 


m-g 


1281 


1264 


m- tl 


568 


260 


m h 


682 


586 


™~ bl 


567 


555 


m^o 
*i 

m x° 2 


87 
325 


84 
415 


m~± 

Xl 


336 


433 


321 


413 


m h 


125 


125 


Br(t2 -> h + h) 





47 


Br(ti -+xi+h+t) 


52 





Br (ft -> 66) 


61 


61 



TABLE I: Two benchmark scenarios. The weak-scale mass 
spectrum and decay branching ratios are calculated using 
SUSY-HIT 0, with m h = 125 GeV set as an input. 

At the LHC, the resonance search of the SM-likc Higgs 
boson via multiple 6-jets has three main backgrounds. 
The first one is the SM background, mainly ttbb. The 
tt+ < 2j is potentially important, but the 46-tagging re- 
quirement can remove much of this background; this is 
discussed further below. Other SM processes involving 
> 46 jets in the final state are less important for our 
purposes; e.g. we have checked that contributions from 
QCD 6666 can be efficiently removed with a $ T and Ht 
requirement, and that ZZ — > 46 makes a negligible con- 
tribution to our search. The second type of background 
is SUSY events containing multiple 6-jets, but with no h 
produced. The third background is combinatorial, aris- 
ing from events where a SM-like Higgs boson decaying 
to 6 quarks is present in the cascade, but where the 6- 
tagged jets paired to reconstruct the Higgs are chosen 
incorrectly. We stress, however, that the latter two types 
are background for the Higgs reconstruction only - in fact 
these 'backgrounds' constitute a strong signal for SUSY. 

The cross sections for all processes here include a K 
factor; Prospino2 |l5| determined this for the SUSY pro- 
cesses. The K factor for pure QCD ttbb and tt+jets is 



2.3, while that for contributions to ttbb from tiZ and tth 
is 1.6 The next-to- leading-order (NLO) cross sec- 

tion of stop pair production at the 8 TeV and 14 TeV 
LHC is shown in Fig. [5] 




FIG. 2: The NLO cross section of stop pair production at the 
8 TeV and 14 TeV LHC generated by Prospino2 [H]. 

Our analysis framework is as follows, 
events with MadGraph.5 / MadEvent 



We generate 
18| , and perform 
showering and hadronization with Pythia 8.1 [19(. For 
the SM backgrounds, we use the Pythia-PGS package in 
MadGraph to do showering and MLM matching. We do 
not employ any dedicated detector simulation in our anal- 
ysis, but we do place cuts on the hadron-level kinematics 
to mimic the response of a generic detector. We subject 
the 4-momentum of each visible final state particle to a 
Gaussian smearing, as implemented in Delphes 2 . [io| . 
After smearing, all visible final state particles are re- 
quired to have px > 0.9 GeV, and \ri\ < 5.0 to fall within 
the calorimeter acceptance. The remaining 'tracks' are 
identified using generator-level information. 

Jet reconstruction is performed next, using 
We use the inclusive anti-fc^ al- 
0.5 and px > 20 GeV, performing 
This sometimes results 



Fast Jet 3.0 [21 1 



gorithm with R 
the clustering using all tracks, 
in 'jets' which are actually isolated leptons and photons. 
Electrons and photons are identified as isolated if the 
scalar sum of charged track pr in a cone of R = 0.2 
does not exceed 10% of the electron or photon pr; 
we further require that the electrons have pr > 20 
GeV and |?7| < 2.47, while isolated muons are instead 
required to have pr > 10 GeV and |7y| < 2.4, with the 
scalar px sum of charged tracks within R = 0.2 of the 
muon less than 1.8 GeV [2(|. Once isolated leptons 
and photons are identified, they are removed from the 
collection of clustered objects. The remaining jets are 
then flavor-tagged with generator-level information. 

B tagging is done by first determining the parton 
with highest pr within the jet cone, and considering 
this to be the 'Monte-Carlo-true' flavor of the jet. True 
6-jets are tagged with an efficiency parameterized by 
e = 0.6tanh(p T /36.0)(1.02 - 0.02|r?|) Q, and true c- 
jets and hght jets are mistagged as 6-jets at flat 10% and 
1% rates [23j, [24|, respectively. Jets are only considered 
for 6-tagging if they fall within the tracker acceptance of 
|t?| < 2.5. Finally, each final-state object is assigned an 
ancestor heavy particle according to the generator-level 



3 



decay history; jet ancestry is determined by the ancestry 
of the hardest parton in the jet cone. 

After the final state objects have been reconstructed, 
events are subjected to selection cuts. Our analysis im- 
plements no triggering step, but it is useful to consider 
this issue briefly The current searches at ATLAS [25| and 
CMS j2(| for tih (which has similar final states and back- 
grounds to our scenarios) use a low-level trigger based on 
isolated leptons. However, such triggers do not have a 
high acceptance for our benchmark points. The focus of 
our strategy is on heavy flavor; it is possible to perform 
6-tagging at the trigger level, so one might consider a 
trigger on multiple 6-tags. Actually, at 8 TeV b-tagging 
is already used at the trigger level, although more as a 
last resort j27|- We would highly suggest explicit multi- 
object triggers with b-tagging at 14 TeV. Alternatively, 
simpler Ht and f T triggers help capture the multiple- 
jet and SUSY nature of our events without the need for 
intensive online computation. Ht > 500 GeV has a high 
efficiency for our benchmark points at 14 TeV; $ T has 
a lower efficiency, but a modest requirement can help re- 
duce much of the QCD background. As a comparison, 
at 8 TeV hadronic p T sum (300 GeV) + MET require- 
ment (100 GeV) is assumed in the CMS experiment fijj ]. 
But, pileup potentially may make these simple triggers 
less robust. 

For the first benchmark point, we apply the following 
cuts at 14 TeV: (1) at least 6 jets with \r]\ < 2.8, (2) at 
least 4 6 tagged jets, at least one with pt > 30 GeV, (3) 
$ T > 150 GeV, and (4) H T > 500 GeV. For the sec- 
ond benchmark point, we amend the last two cuts to (3) 
$ T > 120 GeV, and (4) H T > 650 GeV. In addition to 
the QCD multiple-jets, the 46-tag requirement can also 
efficiently remove the ttbb background. The two non-top 
b quarks in ttbb events are mainly generated by gluon- 
splitting. This induces these events to fail the 46-tag 
cut for two reasons. The first one is at the level of MC 
truth. The b quarks from gluon splitting tend to be softer 
than others, and so recoil more dramatically during par- 
ton showering. This makes the hard-process 6-quarks less 
aligned with the resulting jet, leading to a larger failure 
rate in parton-jet flavor matching. In addition, 6-quarks 
from gluon splitting tend to be more collimatcd, the col- 
limation increasing with the pr of the bb pair; this can 
be seen at parton level in Figure |3] This effect causes the 
resulting jets to overlap and be reconstructed and tagged 
as a single b jet. These effects are implicitly indicated in 
the cut flows of Table HH 

Next, we reconstruct the Higgs resonance from the 6- 
jets in the remaining events, with the Higgs mass strictly 
assumed to be unknown. At this stage, we encounter a 
combinatorial problem. Given that this is a dedicated 
resonance search, we define the signal sample as the cor- 
rectly reconstructed Higgs bosons or the correctly se- 
lected 6-jet pairs. If a wrong 6-jet pair is selected as the 
Higgs candidate for a Higgs event, even if its invariant 




P T (GeV) 

FIG. 3: pt vs. AR for non-top bb quark pairs in ttbb events. 

mass falls into the yet-to-be-known Higgs mass bins, we 
will count it as "combinatorial background" instead of 
a particle signal. This is different from the definition of 
the signal sample in the studies on the SM Higgs search 
via the tth process 23|, 25| and the SUSY non-standard 
Higgs search via the processes associated with additional 
6-jets [2§|, where the Higgs event with a wrongly selected 
6-jet pair as the Higgs candidate, if the 6-jet pair has a 
"correct" invariant mass, was counted as signal instead 
of background. 

The local significance of a resonance is simplest to in- 
terpret when only one dijet per event lies in the local 
mass window Amh, but this method effectively may re- 
sult in the loss of some correct 6-pairs. Multiple pairs 
can be chosen per event as long as their invariant masses 
differ by more than Am/j. To achieve this, we rank the 
jet pairs according to various variables, always keeping 
the first-ranked pair. The second-ranked pair is included 
if |toi — m,2\ > Amj,. and the n tft -ranked pair is included 
if \mi — m„| > Ato/i with i running over all 6-jet pairs 
included before. Am^ is mainly controlled by the reso- 
lution for jet reconstruction at the LHC - in our analy- 
sis, Amh = 40 GeV is assumed. Though an event may 
contribute more than one pair of 6-jets for the resonance 
reconstruction, within the mass window all 6-jet pairs are 
uncorrelatcd. 
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p rank 



FIG. 4: pt vs. # e ff ranking for correct (left) and wrong (right) 
6-jet pairs falling inside the mass window [100, 140] GeV in 
signal events, Case I at 14 TeV. 

To rank the 6-jet pairs in each event, the kinematics of 
jet pairs and jet superstructure [3(| are employed. We use 
the pr of each 6-jet pair, as well as the 'pull angle' of the 
pair. The 'pull' of a jet is the vector in the y — if plane de- 
fined by t = X)ieiet M^^i where n = (Ay ll Aifi) is the 
displacement of the i jet component from the jet axis. 
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t is a measure of the hadronic energy gradient within the 
jet, and carries information about how the jet's ancestral 
parton hadronized. In particular, a pair of jets originat- 
ing from quarks pair-produced by a color singlet tend to 
hadronizc together, so the jet pulls tend to point toward 
each other. For a pair of selected 6-jets with transverse 
momenta p\^ ' b2 , we define an effective "pull angle" 6 e s 
by using the pull angles 9 bl ' b2 which the pulls of the two 
&-jets make with the chord joining the two jets in the 

y - tp plane: 9 cS = ((^/a^)) 2 + (^/a^)) 2 ) ^ 

Here cr(p^ ) = ap^ + b reflects the jet px dependence of 
the standard deviation of the pull angle for the two &-jets 
produced from Higgs decay, as noticed in [3(J. For the 
kinematic regime that we are considering, a = — 1 TeV -1 
and b = 1.5 are assumed. 

If the 6-jet pairs in each event arc ranked according to 
these variables, noticeable differences emerge. Figure 0] 
shows the pt vs. 9 c g ranks of correct and wrong pairs 
in Case I at the 14 TeV LHC. The distribution of pairs 
in the ranking plane can be used to improve the Higgs 
search sensitivity. 



y/s = 14 TeV 
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Case I 


Case II 


Events 


5.2 x 10' 


8.2 x 10 s 


26176 
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Cut 1 
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474234 
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Cut 2 


88700 


12077 


961 


824 


Cut 3 


51 / 79 


442 / 796 


567 


436 


Cut 4 


29 / 23 


351 / 366 


547 


389 


Choice A 


20 / 11 


5+157 / 4+157 


99+215 


65+135 


Choice B 


20 / 12 


4+166 / 4+159 


91+219 


86+108 


Choice C 


13 / 13 


5+104 / 3+104 


78+147 


65+62 



TABLE II: Cut flows for the benchmark point at 14 TeV. The 
cut flow of some SM events is labeled as 'Case I / Case II', 
where different cuts are used at the same stage, tibb includes 
QCD, tiZ and tih. Choice A, B, C rows correspond to the 6 c h , 
Pt and pr — d a s plane pair selection methods, respectively. 
Bins labeled as (N + M) show (Higgs + wrong) 6-jet pairs. 




m bb (GeV) 



FIG. 5: B-jet pair invariant mass, Case I at 14 TeV. 

Complete cut flows of both benchmark points at 14 
TeV are presented in table [TTJ The number of events have 
to the cross section of each process. A first observation is 
a relatively large sensitivity to natural SUSY in our cho- 
sen benchmark points. From Table |H] the sensitivities 
in both cases are > 20a for 40fb _1 data at the 14 TeV 




FIG. 6: Signal significance using various methods of choosing 
pairs, Case I at 14 TeV, with 40 fb -1 of data. 




m bb (GeV) 



FIG. 7: B-jet pair invariant mass, Case II at 14 TeV. 




FIG. 8: Signal significance using various methods of choosing 
pairs, Case II at 14 TeV, with 40 fb" 1 of data. 

LHC. Again, note that the contribution labeled "SUSY: 
wrong pair" in these figures is a 'background' only for the 
Higgs search, but indicates a large signal for the discov- 
ery of SUSY itself. To combine the pr and 6 c ff ranking 
strategies, we require that the selected pairs fall within 
the triangular region where rank(pr) + rank (# c ff) < 5. 
This strategy is used for the invariant mass plots of the 
selected 6-jet pairs in Figure [S] and Figure [71 for the SM- 
likc Higgs candidates. The local significances for different 
strategies of the resonance reconstruction are shown in 
Figure [n]and Figure [HI For the best option, the correctly 
paired Higgs 6-jets give a local significance of > 5cr in 
both cases for 40 fb -1 of data, with the local mass win- 
dow taken to be 100 GcV to 140 GeV. Though ranker) 
+ rank (9 e g) < 5 does not increase the values of S/V~B, 
it improves the S/B to be above 0.3 and 0.4 for the two 
cases, respectively. This can decrease the impact of sys- 
tematic uncertainties. Each of these strategies is superior 
to naively selecting the b pair with invariant mass closest 
to 125 GeV; for both benchmark points, this naive choice 
gives a sensitivity of ~ 4cr in Case I and ~ 3.5<r in Case 
II, with S/B ~ 0.17 in both cases. 
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Our ^-tagging effi- 
If the b-tagging ef- 



These results may be improved, 
ciency is rather conservative 31 1. 
ficiency is assumed to be 0.7, similar to the approach 
in 0, HH , then the significances for various strategies are 
expected to be universally increased by a factor ~ 1.5. 
Also, additional kinematic variables, such as "twist an- 
gle" [33l ] , may improve on the performance of those used 
in our analysis. 

Although the foregoing discussion is confined within a 
special theory beyond the SM - natural SUSY - its im- 
pact is fairly profound. It provides an independent way 
to discover the SM-like Higgs boson and to understand 
the mass origin of the SM fermions. In turn, searching for 
the Higgs boson using this strategy provides a direct way 
to test natural SUSY. In addition, the potential applica- 
ble scope of this strategy is very broad. The first class 
of examples are some non-SUSY theories with naturally 
EWSB, such as little Higgs models [34j . composite Higgs 
models 351, and Randall-Sundrum models 13611. where 



[6] 



fcrmionic top partners t 1 with a mass below 1 TeV are 
typically predicted [37j . given a 125 GeV SM-like Higgs 
boson, and the Higgs boson can be produced via t' — > ht, 
a way similar to the second mechanism in natural SUSY. 
Actually, multiple 6-tagging has been noticed to be useful 
in this context [38j], while the Higgs reconstruction was 
thought to be a big problem due to the combinatorial 
background [39|. Another class of examples are the SM 
Higgs search via the tih production and the non-standard 
Higgs (H, A) search via the bbH or bbA production in 
the MSSM. Though the involved kinematics are differ- 
ent, given h,H,A — > bb, these analyses are expected to 
share the feature of suppressed 4j and tt, bb + 2j back- 
ground, leaving the continuum ttbb or bbbb as the main 
one. Then the Higgs resonance can be reconstructed with 
color-flow or kinematic variables. We leave the consider- 
ation of these interesting possibilities to future work. 
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